The kinetic and interaction energies of a trapped Bose gas: Beyond the mean field 
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The kinetic and interaction energies of a three-dimensional (3D) dilute ground-state Bose gas 
confined in a trap are calculated beyond a mean-field treatment. They are found to depend on 
the pairwise interaction trough two characteristic lengths: the first, a, is the well-known scattering 
length and the second, b, is related to the latter by 6 = a — \da/d\ with A being the coupling 
constant. Numerical estimations show that the pairwise interaction energy of a dilute gas of alkali 
atoms in a trap is negative (in spite of the positive scattering length); its absolute value is found by 
about the order of magnitude larger than that of the mean-field interaction energy that corresponds 
to the last term in the Gross-Pitaevskii (GP) functional. 
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It is known that the total energy i? of a dilute 3D 
Bose gas confined in a trap can be represented at T ^ Tc 
(Tc is the temperature of the Bose- Einstein condensation) 
as the functional of the order parameter (j> = (f'i^) = 
('0(r)) (here ?/'(r) is the Bose field operator) 



:^|V^P + K,t(r)|^P + 
2 m 



\<t>\' 



(1) 



called the Gross-Pitaevskii functional. In Eq. (|^) V^xt (r) 
denotes an external trapping field with the characteristic 
length L, typical for the spatial variation of this field [for 

a harmonic trap, L ~ aho — (mujho) ], and |a| <C i 
stands for the scattering length in a gas considered. From 
the definition it follows that the order parameter satis- 
fies the normalization condition J c?'^r|(/)p = N, where the 
number of bosons Nq in the Bose-Einstein condensate is 
replaced by their total number N due to a small conden- 
sate depletion. The stationary solution corresponding to 
the minimum of the functional (|l|) obeys the relations 

6{E - fiN)/S(l)*{r) ^ 5{E - fiN)/S(j){r) = 0, (2) 



which give the GP equation 



^0 : 
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+ Kxt(r)<^-l- 1<^|'0. (3) 



The term fiN appears in Eq. (|2|) due to the normaliza- 
tion condition, and the Lagrange multiplier is nothing 
but the chemical potential. In this paper we address the 
problem of interpretation of different terms in the GP 
functional; namely, the kinetic ii^kin md pairwise inter- 
action Eint energies are under investigation below. Note 
that this is of actual interest because analysis of the ex- 
periments with magnetically trapped Bose gases is car- 
ried out usually in terms of Skin and Eint j^i^- More- 
over, the problem of that interpretation is rather ambigu- 
ous due to just opposite points of view on this subject 
found in the literature. Indeed, in the paper |^ the third 
term of the GP functional is treated as the energy of the 



pairwise boson interaction. While according to another 
point of view Q , in the homogeneous case 14xt = this 
term makes contribution only to the kinetic energy if the 
boson-boson interaction potential V{r) is of the hard- 
sphere form: V{r) = +oo at r < a, and V{r) = at 
r ^ a. 

This problem is usually considered to be trivial 
Indeed, in the mean-field interpretation ||^, the terms in 
Eq. (j^) are associated with the kinetic energy, the energy 
of interaction with the external field, and the pairwise 
interaction energy, respectively. However, in the partic- 
ular case T4xt — 0, this interpretation contradicts both 
the results of Lieb and Yngvason Q and ours To 
have an understanding of the situation, let us calculate 
the kinetic and interaction energies of a trapped Bose gas 
on a solid theoretical basis beyond the mean-field inter- 
pretation of the terms of the GP functional. This can 
be realized with the help of the well-known variational 
theorem for the grand canonical potential at zero tem- 
perature Sn = 6{E - fj.N) = {6{H - nN)), where (• • •) 
stands for the statistical average over the grand canonical 
ensemble. To calculate the pairwise interaction energy 
Eint = ^(ki ~ i"jl)/2), it is useful to introduce the 

coupling constant A [V{r) XV (r), and A = 1 in final 
formulas] . Then the thermodynamic variational theorem 
leads to Eint — d{E — fiN)/dX, which, in conjunction 
with Eqs. (0) and (|), yields § 



Eu 



2ntf da 
m dX 



(4) 



where (f) is the stationary order parameter obeying the 
GP equation (^). Being the characteristic of the two- 
body problem, the derivative da/dX is obtained from the 
other variational theorem proved in the papers 1^ for the 
short-range pairwise interaction potentials that go to zero 
at r ^ oo as V{r) — * l/r™ (m > 3), or even faster. The 
particular variant of this theorem allows for connecting 
the infinitesimal change of the scattering length 5a with 
that of the potential 6V{r): 
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d'V[(^(")(r)]2^T/(r), 



(5) 



where (^'''^(r) is the s-wave solution of the two-body 
Schrodinger equation in the center-of-mass system 



(6) 



Equation (^) corresponds to the scattering state with the 
momentum p — 0, and its solution ip^^^ (r) is chosen to be 
real due to the real boundary condition ip^^\r) l~a/r 
at r oo. This boundary condition and Eq. (|^) yield 
the following relation: 



rf3rF(r)v3W(r) 



Equations (||) and (g) lead to 

where we put by definition 

6^ i_ /■ dV|V(^W(r)|'. 



(7) 
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From Eq. (|^) it follows that 6 is a positive quantity and 
can be considered as a new characteristic length, which is 
not expressed in terms of a and depends on a particular 
shape of the interaction potential V{r). For example, for 
the hard spheres [V{r) — +oo at r < a, and V{r) = at 
r ^ a] we have 6 = a ||^. While for V{r) close to zero, 
in the weak-coupling regime we obtain b < \a\ [§. 
Thus, Eq. (0), taken together with Eq. (0), gives 



Eint = (a - b) 

m 



d^r\(l3\^ 



(10) 



In the same manner, using the thermodynamic varia- 
tional theorem with respect to AVcxt(r), one obtains for 
the energy of interaction with the external field 



Eext — 



derivation of the GP functional jg] beyond the mean-field 
approximation. Worth mentioning that this approxima- 
tion has recently been demonstrated |^ to face significant 
problems in the strong-coupling regime The point 
is that the mean-field approach yields the correct result 
for the total energy of a Bose gas of strongly interacting 
particles (strictly speaking, after removing the ultraviolet 
divergence [||). However, it fails to give the appropriate 
picture of the short-range boson correlations in this case, 
which leads to the ultraviolet divergence and incorrect 
values for the kinetic and interaction energies [^. We 
stress that this only concerns the strong-coupling case. 
There is no problem with the mean-field theory in the 
weak-coupling regime, when b/\a\ <C 1, and, hence, one 
can approximately put 6 = in Eqs. (^0|) and Note 
that in the uniform case V^xt — the GP equation (^) 
gives (j> ~ const = \/n (here n = N/V is the density of 
bosons), and Eqs. ( p^ ) and ( p^ are reduced to the results 

Eint/N = 27rfi^n(a - b)/m, Ey^i^/N = 2T:H'^nb/m, 

derived earlier by the present authors 

Mathematically rigorous, the derivation 1^ of the GP 
functional does not give any information about nature of 
the third term in Eq. ( fi\). This nature becomes, of course, 
clear from Eqs. (pO[)-(|12D, but the method of obtaining 
them is also rather formal. For a more deep insight, let us 
find out what approximations for the reduced density ma- 
trices result in Eqs. ( |l0[ ) and (^l]). Introducing the one- 
body (1-matrix) Fi(r, r') = ('i/;^(r)'(/;(r')) and two-body 
(2-matrix) i^2(ri, ra; r'^, r'2) = {i^Hn)^;^ {r2Mr'2Mr[)) 
density matrices, one writes 

^^ext = y"dVFcxt(r)Fi(r,r), (13) 

Eint^l j d\d\'V{\v-v'\)F2{v,v'-vy). (14) 

The 1-matrix can be expanded in the orthonormal set of 
its eigenprojectors as 



|^V;,,(r,)) - j dVK=xt(r)|0r (11) 



7Vo0S(r)0o(r') 



n.^:(r)0,(r'). (15) 



In turn, the kinetic energy can be calculated by using 
Eqs. (|l|), 1^ and (|ll|) with the resuh 



2m/ 



2TTn%. 



(12) 



Thus, the thermodynamic variational theorem and mean- 
field treatment [the latter leads to Eqs. (p^)-(p^) but with 
6 = 0] yield in general different results for the kinetic and 
pairwise interaction energies of a trapped Bose gas. Here 
one should not be confused by the fact that the GP equa- 
tion is usually considered to be a product of the mean- 
field approach. Though it has originally been derived 
by means of this approach [Q, there exists the rigorous 



Here (/)o(r) — 0(r)/-yiVo = {tp{r))/^/NQ is the eigen- 
function corresponding to the macroscopic eigenvalue Nq 
[/dVFi(r',r)0o(r') = ^o0o(r), / d^r |0o(r)P = 1], see, 
e.g., Ref. By definition, the one-body matrix is nor- 
malized as / d^r Fi{r, r) = N; therefore, in the case of a 
small condensate depletion {N ~ No)/N <C 1 one can put 

Fi(r,r')=^(^nr))(V;(r')), 

which, taken together with Eq. (^3|) , results immediately 
in Eq. (|ll]) . Note that the calculation of the second term 
of the kinetic energy ( p^ ) which involves the 1-matrix 
requires knowledge of the other eigenvalues and eigen- 
functions of i^i(r,r'). However, this calculation is not 
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really needed because, when treating a trapped Bose gas 
on the basis of the microscopical information concerning 
the one-body and two-body matrices, we are able to re- 
strict ourselves only to derivation of Eqs. (^0|) and (pi]). 
As to Eq. (p^), it can be determined from the relation 
E = E^n + E,^t + Eint with the help of Eqs. (|l|), 
and (fill) . Let us repeat one more that the second term in 
Eq. ( |l2|) is exactly due to the second term in Eq. (15|). For 
this reason, the approximation for the momentum distri- 
bution Up ~ |(/>(p)p (where (f>{p) is the Fourier transform 
of the order par ameter) is able to produce only the first 
term in Eq. @. 

Explicit expressions for a number of eigenvalues and 
eigcnfunctions of the two-body matrix have been ob- 
tained for a homogeneous system of bosons in the pa- 
per and were then successfully employed in the 
papers |^. In particular, the anomalous two-operator 
(quasi-) average {■i(:{ri)iip{r2)) / \/ No{No — 1) was shown 
to be the eigenfunction (pair wave function) of the 
2-matrix, corresponding to the maximum macroscopic 
eigenvalue No{No — 1) ~ Nq, which is nothing else 
but the number of pairs of bosons in the Bose-Einstein 
condensate. This result is also valid for a non-uniform 
(trapped) Bose system. In the homogeneous case, this 
avera ge, contrary to the situation for the Fermi sys- 
tems |ll[| , describes the scattering (not bound) state of 
a pair of the condensed bosons. From this interpreta- 
tion it follows that at small separations the pair wave 
function (^(ri)'0(r2))/v/A^o(A^o - 1) should be propor- 
tional to the two-body wave function (p^"-* (r) obeying the 
Shrodinger equation On the other hand, the Bogoli- 
ubov principle of the correlation weakening [|ll[ yields 

(^(ri)^(r2)) ('0(ri))('0(r2)) when |ri - oo. 
From the physical point of view, the latter limit implies 
|ri — r2| » Icoh, where the length of coherence |^ /coh = 
l/y/Sirna (here n ~ |0(ri)p ~ |(/)(r2)P is the local den- 
sity) is assumed to obey the inequalities |a| < Zcoh <C L. 
Therefore, we arrive at the following approximation: 

(V'(ri)V'(r2)) ^ 0(ri)0(r2)^(")(|ri - r2|). (16) 

Remind that (p'^^^{r) is the solution of Eq. (^) with the 
asymptotics (p^-^'){r) — > 1 — a/r for r — s- oo. We stress 
that the approximation (|l^) works well at |ri — r2| < 
Icoh- By analogy with the 1-matrix, the small condensate 
depletion leads to the representation 

F2(ri,r2;r;,r'2) ~ {^p\r^)^/J^ (r^)) {^lj{r'2)4^{r[)) , (17) 

which, with the help of Eqs. ( [T^ ) and (p^), yields 

Eint d\d\' V{\v - r'|)[(^W(|r - r'l)]^ 



In order to estimate the ratio b/a typical for the alkali 
atoms, we make use of the model potential 



x|0(r)n0(r')| = 



(18) 



Finally, using the inequality L ^ a valid in the GP ap- 
proximation and the theorem (^) at A = 1, we arrive at 
Eq. (0). 



V{r) 



-oo. 



r > ro- 



(19) 



It is usually considered (see, e.g., ||T^) to be relevant in 
the situation of interest. This model interaction leads jl^ 
to the scattering length 



a/r, = r(3/4)J_i/4(a;o)/[2r(5/4)Ji/4(a;o)]. 



(20) 



Here xq = r1/{2r1); J^{x) and T{z) denote the Bessel 
function and the Euler gamma- function, respectively. 
Further calculating da/dX with Eq. ( pO| ) Ir, X^^^r, 
in Eq. ([19|)] and, then, using Eqs. (||) and ( po|) at A = 1, 
we arrive at 

b/a = 3/4 + 1/ [7rV2 Ji/4(xo) J_i/4(xo)] • (21) 



For the potential (^^, the parameter xq is connected 
with the number of the bound states for the two-body 
problem: larger the number, larger xq- It is well-known 
that this number is equal to that of the nodes of the 
wave function (p^°\r) obeying Eq. (|[). The pairwise in- 
teraction of alkali atoms is characterized by large amount 
of possible bound pair states, which, however, kinetically 
inaccessible in the trapped systems [Q. So, we are inter- 
ested in the regime when xq ^ 1. According to Eq. (|2l|), 
b/\a\ » 1 at xo » 1. The same situation {b/\a\ » 1) 
takes place also at any xq provided the denominator in 
Eq. ( [20| ) is close to zero (as in the case of the Feshbach 
resonance). The large value of b/\a\ means that the sec- 
ond term in Eq. (|l2|) is positive and much larger than 
the absolute value of the third term in the GP functional, 



while the pairwise interaction energy ( 10 ) is negative and 
close, in absolute value, to the second term of Eq. (p^). 
In other words, for a homogeneous dilute Bose gas of al- 
kali atoms [when the first term in Eq. ( p^ equals to zero] 
we can expect i?kin — \Eint\ ^ \E\. In order to obtain 
numerical estimations for sodium atoms, we utilize the 
data of Ref. ||lj| r, = 88. lap and a = 52ao (oq is the 
Bohr radius). Note that, at any r, and a, solutions of 
Eq. ( po[ ) for the core radius rp form an infinite sequence, 
approaching zero; each of the solutions for ro corresponds 
to a certain number A^nodc of the nodes of ip'^^\r). Nu- 
merical results for b/a versus A^node are shown in Fig. 
|l|(a). A typical solution of the Schrodinger equation (g) 
with the potential ( p^ corresponding to A^nodc = 12 [and, 
hence, = 10.03ao] is given in Fig. |l|(b). Now, the Na- 
Na interaction is expected [^ij to provide more than 15 
possible bound states, and, thus, we have the estimation 
b/a > 20 [see Fig. |](a)]. We remark that the visible 
linear dependence of b/a on the number of the nodes is 
a specific feature of the attractive interatomic potential 
proportional to 1/r^ but not a general relation between 
b and a. We also stress that the pair distribution func- 
tion of sodium atoms given in Fig. 0(c) is equal to zero 
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FIG. 1. (a) The ratio hja versus the number A^nodc of the nodes of (^'''•'(r) [(^'''•'(r) is the s-wave solution of the 
Schrodinger equation with V(r) given by Eq. (|l9|)] for the given values of the scattering length a = 52ao and the 
model parameter Vc = 88.1ao. This ratio is equal to E]^i-a/E for a homogeneous dilute Bose gas of sodium atoms, 
(b) The wave function (^^'"^ (r) for the model potential ( [l9| ) with N-^odc = 12. (c) The pair distribution function 
g(|ri - r2|) = F2(ri, r2; ri, r2]/[|(?!)(ri) p |(?!)(r2) ~ [tp'^°\\vi - V2\)f [see Eqs. (|l|) and Q] in a trapped sodium atoms, 
interacting via the potential ([19|) with A'^„odc ~ 5. Intensive oscillations occur when r < a — 52ao. 



at r = 0, contrary to the expectations |^ inspired by the 
mean- field treatment. 

In conclusion, the kinetic and interaction energies (p^- 
12|) of a 3D dilute ground-state Bose gas confined in a 
trap have been derived beyond the mean-field interpreta- 
tion. The pairwise interaction energy i?i„t is found to be 
controlled by the two characteristic lengths: in addition 
to the well-known scattering length a [Eq. (0)], the ex- 
pression ( p^ ) involves the positive parameter b [Eq. (^]. 
Whereas the kinetic energy (|l^) depends on the pairwise 
interaction only through the characteristic length h. The 
derived estimations suggest that in the experimentally 
interesting case of a dilute Bose gas of alkali atoms the 
part of the kinetic energy coming from the pairwise in- 
teraction is much larger than the absolute value of the 
third term in the GP functional (|l|) . 

Note that experimental observation of i?kin and i^int 
is quite possible provided the diagonal element of the 2- 
matrix F2(r, r'; r, r') (the familiar pair distribution func- 
tion) and the momentum occupation number Up — 
(ajjOp) are experimentally determined. Indeed, the pair- 
wise interaction energy is then calculated by means of 
Eq. (|l|) with the pairwise potential V{r) that can be 
found numerically (see, e.g., Ref. [|T^). The kinetic en- 
ergy per particle can be represented as 



(27r)3 2m ~n 



(22) 



The quantities F2 and Up are the well-known subjects 
of experimental investigations in the condensed matter 
physics, and there is technique of measuring these quan- 
tities (see, e.g., discussion in Ref. [|l^). However, in 
the situation of interest, the trouble point is that, in or- 
der to calculate -Bkm, one needs to know Up beyond the 
phonon regime for p > H/a because the contribution of 
the phonon region is small in the integral in Eq. (12? 



Moreover, measurements with sufficiently small step in 
the momentum direction are of importance here to in- 
tegrate properly in Eq. (|22). The same is related to 



Eq. (p^, where the short-range behaviour |r — r'| < a 
of the pair distribution function F2(r, r'; r, r') is of im- 
portance. For this reason, accurate data are needed for 
the static structure factor at p > H/a. Moreover, one 
again needs to make fine step-by-step measurements. To 
the best of our knowledge, so far there are only several 
experimental points for the static structure factor [^. 
The situation concerning Up is now even worse. Thus, 
additional experiments are needed to make final conclu- 
sions about the failure of the mean-field approach. 

It is also worth noting that, in experiments on the ex- 
pansion of a condensate, one observes a kinetic energy of 
the freely expanding condensate, which is equal to a sum 
of the pairwise potential (10) and kinetic ( [l2| ) energies of 
the trapped one, in effect proportional to the scattering 
length. 
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